The Influence of Snow Cover Variability and Tundra Lakes on Passive Microwave Remote Sensing of Late Winter Snow Water Equivalent in the Hudson Bay Lowlands

Current North American operational satellite passive microwave snow water equivalent (SWE) retrieval algorithms consistently underestimate SWE levels for tundra environments when compared to four years of regional snow surveys conducted in the Northwest Territories and northern Manitoba, Canada. Almost all contemporary SWE algorithms are based on the brightness temperature difference between the 37GHz and 19GHz frequencies found onboard both past and present spaceborne sensors. This underestimation is likely a result of the distribution and deposition of the tundra snow, coupled with the influence of tundra lakes on brightness temperatures at the 19GHz frequency. To better our understanding concerning the underestimation of passive microwave SWE retrievals on the tundra, Environment Canada collected in situ measurements of SWE, snow depth, and density at 87 sites within a 25km by 25km study domain located near Churchill, Manitoba in March 2006. Coincident multi-scale passive microwave airborne (70m & 500m resolution) and spaceborne (regridded to 12.5km & 25km resolution depending on frequency) data were measured at 6.9GHz, 19GHz, 37GHz and 89 GHz frequencies during the same time period. The snow survey data highlighted small-scale localized patterns of snow distribution and deposition on the tundra that likely influences current SWE underestimation. Snow from the open tundra plains is re-distributed by wind into small-scale vegetated features and micro-topographic depressions such as narrow creekbeds, lake edge willows, small stands of coniferous trees and polygonal wedge depressions. The very large amounts of snow deposited in these spatially-constrained features has little influence on the microwave emission measured by large-scale passive microwave spaceborne sensors and is therefore unaccounted for in current methods of satellite SWE estimation. The analysis of the passive microwave airborne data revealed that brightness temperatures at the 19GHz were much lower over some tundra lakes, effectively lowering SWE at the satellite scale by reducing the 37-19GHz brightness temperature difference used to estimate SWE. The unique emission properties of lakes in the wide open expanse of the tundra plains, coupled with an insensitivity to the large amounts of SWE deposited in small-scale features provides an explanation for current passive microwave underestimation of SWE in the tundra environment

Snow stratigraphic heterogeneity within ground-based passive microwave radiometer footprints: implications for emission modeling

Two-dimensional measurements of snowpack properties (stratigraphic layering, density, grain size and temperature) were used as inputs to the multi-layer Helsinki University of Technology (HUT) microwave emission model at a centimeter-scale horizontal resolution, across a 4.5 m transect of ground-based passive microwave radiometer footprints near Churchill, Manitoba, Canada. Snowpack stratigraphy was complex (between six and eight layers) with only three layers extending continuously throughout the length of the transect. Distributions of one-dimensional simulations, accurately representing complex stratigraphic layering, were evaluated using measured brightness temperatures. Large biases (36 to 68 K) between simulated and measured brightness temperatures were minimized (-0.5 to 0.6 K), within measurement accuracy, through application of grain scaling factors (2.6 to 5.3) at different combinations of frequencies, polarizations and model extinction coefficients. Grain scaling factors compensated for uncertainty relating optical SSA to HUT effective grain size inputs and quantified relative differences in scattering and absorption properties of various extinction coefficients. The HUT model required accurate representation of ice lenses, particularly at horizontal polarization, and large grain scaling factors highlighted the need to consider microstructure beyond the size of individual grains. As variability of extinction coefficients was strongly influenced by the proportion of large (hoar) grains in a vertical profile, it is important to consider simulations from distributions of one-dimensional profiles rather than single profiles, especially in sub-Arctic snowpacks where stratigraphic variability can be high. Model sensitivity experiments suggested the level of error in field measurements and the new methodological framework used to apply them in a snow emission model were satisfactory. Layer amalgamation showed a three-layer representation of snowpack stratigraphy reduced the bias of a one-layer representation by about 50%

Evaluation of Operation IceBridge quick-look snow depth estimates on sea ice

We evaluate Operation IceBridge (OIB) ‘quick-look’ (QL) snow depth on sea ice retrievals using in situ measurements taken over immobile first-year ice (FYI) and multi-year ice (MYI) during March of 2014. Good agreement was found over undeformed FYI (-4.5 cm mean bias) with reduced agreement over deformed FYI (-6.6 cm mean bias). Over MYI, the mean bias was -5.7 cm but 54% of retrievals were discarded by the OIB retrieval process as compared to only 10% over FYI. Footprint scale analysis revealed a root mean square error (RMSE) of 6.2 cm over undeformed FYI with RMSE of 10.5 cm and 17.5 cm in the more complex deformed FYI and MYI environments. Correlation analysis was used to demonstrate contrasting retrieval uncertainty associated with spatial aggregation and ice surface roughness

Observing Scattering Mechanisms of Bubbled Freshwater Lake Ice Using Polarimetric RADARSAT-2 (C-Band) and UW-Scat (X- and Ku-Bands)

A winter time series of ground-based (X- and Ku-bands) scatterometer and spaceborne synthetic aperture radar (SAR) (C-band) fully polarimetric observations coincident with in situ snow and ice measurements are used to identify the dominant scattering mechanism in bubbled freshwater lake ice in the Hudson Bay Lowlands near Churchill, Manitoba. Scatterometer observations identify two physical sources of backscatter from the ice cover: the snow-ice and ice-water interfaces. Backscatter time series at all frequencies show increases from the ice-water interface prior to the inclusion of tubular bubbles in the ice column based on in situ observations, indicating scattering mechanisms independent of double-bounce scatter. The co-polarized phase difference of interactions at the ice-water interface from both scatterometer and SAR observations is centered at 0° during the time series, also indicating a scattering regime other than double bounce. A Yamaguchi three-component decomposition of the RADARSAT-2 C-band time series is presented, which suggests the dominant scattering mechanism to be single-bounce off the ice-water interface with appreciable surface roughness or preferentially oriented facets, regardless of the presence, absence, or density of tubular bubble inclusions. This paper builds on newly established evidence of single-bounce scattering mechanism for freshwater lake ice and is the first to present a winter time series of ground-based and spaceborne fully polarimetric active microwave observations with polarimetric decompositions for bubbled freshwater lake ice.European Space Agency (ESTEC): 10.13039/501100000844 Natural Sciences and Engineering Research Council of Canada: 10.13039/50110000003

Radio-frequency interference mitigating hyperspectral L-band radiometer

Radio-frequency interference (RFI) can significantly contaminate the measured radiometric signal of current spaceborne L-band passive microwave radiometers. These spaceborne radiometers operate within the protected passive remote sensing and radio-astronomy frequency allocation of 1400–1427 MHz but nonetheless are still subjected to frequent RFI intrusions. We present a unique surface-based and airborne hyperspectral 385 channel, dual polarization, L-band Fourier transform, RFI-detecting radiometer designed with a frequency range from 1400 through  ≈  1550 MHz. The extended frequency range was intended to increase the likelihood of detecting adjacent RFI-free channels to increase the signal, and therefore the thermal resolution, of the radiometer instrument. The external instrument calibration uses three targets (sky, ambient, and warm), and validation from independent stability measurements shows a mean absolute error (MAE) of 1.0 K for ambient and warm targets and 1.5 K for sky. A simple but effective RFI removal method which exploits the large number of frequency channels is also described. This method separates the desired thermal emission from RFI intrusions and was evaluated with synthetic microwave spectra generated using a Monte Carlo approach and validated with surface-based and airborne experimental measurements

Effect of snow microstructure variability on Ku-band radar snow water equivalent retrievals

Spatial variability in snowpack properties negatively impacts our capacity to make direct measurements of snow water equivalent (SWE) using satellites. A comprehensive data set of snow microstructure (94 profiles at 36 sites) and snow layer thickness (9000 vertical profiles across 9 trenches) collected over two winters at Trail Valley Creek, NWT, Canada, were applied in synthetic radiative transfer experiments. This allowed robust assessment of the impact of estimation accuracy of unknown snow microstructural characteristics on the viability of SWE retrievals. Depth hoar layer thickness varied over the shortest horizontal distances, controlled by subnivean vegetation and topography, while variability of total snowpack thickness approximated that of wind slab layers. Mean horizontal correlation lengths of layer thickness were sub-metre for all layers. Depth hoar was consistently ~30% of total depth, and with increasing total depth the proportion of wind slab increased at the expense of the decreasing surface snow layer. Distinct differences were evident between distributions of layer properties; a single median value represented density and specific surface area (SSA) of each layer well. Spatial variability in microstructure of depth hoar layers dominated SWE retrieval errors. A depth hoar SSA estimate of around 7% under the median value was needed to accurately retrieve SWE. In shallow snowpacks <0.6m, depth hoar SSA estimates of ±5-10% around the optimal retrieval SSA allowed SWE retrievals within a tolerance of ±30 mm. Where snowpacks were deeper than ~30cm, accurate values of representative SSA for depth hoar became critical as retrieval errors were exceeded if the median depth hoar SSA was applied

SLAPex Freeze/Thaw 2015: The First Dedicated Soil Freeze/Thaw Airborne Campaign

Soil freezing and thawing is an important process in the terrestrial water, energy, and carbon cycles, marking the change between two very different hydraulic, thermal, and biological regimes. NASA's Soil Moisture Active/Passive (SMAP) mission includes a binary freeze/thaw data product. While there have been ground-based remote sensing field measurements observing soil freeze/thaw at the point scale, and airborne campaigns that observed some frozen soil areas (e.g., BOREAS), the recently-completed SLAPex Freeze/Thaw (F/T) campaign is the first airborne campaign dedicated solely to observing frozen/thawed soil with both passive and active microwave sensors and dedicated ground truth, in order to enable detailed process-level exploration of the remote sensing signatures and in situ soil conditions. SLAPex F/T utilized the Scanning L-band Active/Passive (SLAP) instrument, an airborne simulator of SMAP developed at NASA's Goddard Space Flight Center, and was conducted near Winnipeg, Manitoba, Canada, in October/November, 2015. Future soil moisture missions are also expected to include soil freeze/thaw products, and the loss of the radar on SMAP means that airborne radar-radiometer observations like those that SLAP provides are unique assets for freeze/thaw algorithm development. This paper will present an overview of SLAPex F/T, including descriptions of the site, airborne and ground-based remote sensing, ground truth, as well as preliminary results

Analysis of Soil Freeze/Thaw Signatures During Slapex F/T Campaign

Permanently frozen and seasonally frozen soils occur over a large portion of the Earth's land surface. Changes in the freeze/thaw state of the land surface reflects major changes in thermal and hydraulic properties as well as acting as a "switch" for many ecological processes. In short, soil freeze/thaw state is a fundamental land surface variable in the water and energy cycles, and it connects to the carbon cycle. Surface freeze/thaw state is observable by passive and active microwave sensors. For example, NASA's Soil Moisture Active Passive (SMAP) mission includes a freeze/thaw data product. Such satellite sensing offers routine all-season and all-weather global observations of soil freeze/thaw state with the application of suitable algorithms. We describe early finding from the SLAPex Freeze/Thaw campaign, believed to be the first airborne campaign of its type, focusing on soil freeze/thaw

Meteorological observation from Trail Valley Creek (TVC, Canada) collected by two automatic weather stations during winter 2012-2013

Meteorological observations for every half an hour from Trail Valley Creek (TVC), Canada, collected by two automatic weather stations; one in the upper plateau and the other in the valley, during the winter 2012-2013